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Abstract 

We present a general numerical method for computing guaranteed two-sided bounds 
for principal eigenvalues of symmetric linear elliptic differential operators. The approach 
is based on the Galerkin method, on the method of a priori-a posteriori inequalities, 
and on a complementarity technique. The two-sided bounds are formulated in a general 
Hilbert space setting and as a byproduct we prove an abstract inequality of Friedrichs'- 
Poincare type. The abstract results are then applied to Friedrichs', Poincare, and trace 
inequalities and fully computable two-sided bounds on the optimal constants in these 
inequalities are obtained. Accuracy of the method is illustrated on numerical examples. 
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1 Introduction 

Eigenvalue problems for differential operators have attracted a lot of attention as they have 
many applications. These include the dynamic analysis of mechanical systems [3, 15, 25], 
linear stability of flows in fluid mechanics [24] , and electronic band structure calculations [29] . 
In this paper, we concentrate on guaranteed two-sided bounds of the principal (smallest) 
eigenvalue of symmetric linear elliptic operators. The standard Galerkin method for solution 
of eigenproblems is efficient and its convergence and other properties are well analysed [5, 
6, 11]. It is also well known for providing upper bounds on eigenvalues. However, in many 
applications a reliable lower bound of the smallest eigenvalue is the key piece of information 
and, unfortunately, the Galerkin method cannot provide it. 

The question of lower bounds on the smallest eigenvalue has already been studied for 
several decades. For example see [30], where the second order elliptic eigenvalue problems 
with Dirichlet boundary conditions are considered. Another technique that gives the lower 
bounds not only for the first eigenvalue is the method of intermediate problems. It is based 



*Thc support of I. Scbcstova by the Grant SVV-2013-267316 and by the project MathMAC - University 
center for mathematical modeling, applied analysis and computational mathematics of the Charles University 
in Prague and the support of T. Vejchodsky by RVO 67985840 are gratefully acknowledged. 

t Department of Numerical Mathematics, Faculty of Mathematics and Physics, Charles University in 
Prague, Sokolovska 83, CZ-186 75 Praha 8, Czech Republic (ivasebestovaOseznam. cz). 

■^Institute of Mathematics, Academy of Sciences, Zitna 25, CZ-115 67 Praha 1, Czech Republic 
(ve j chodOmath .cas.cz). 



1 



TWO-SIDED BOUNDS FOR EIGENVALUES 



2 



on finding a base problem and subsequently introducing intermediate problems that give 
lower bounds for eigenvalues of the original problem and at the same time can be resolved 
explicitly, sec for example [8, 9, 39]. A survey of this technique can be found in [16]. 

Nonconforming methods have been used for computing lower bounds on eigenvalues, see 
for example [4, 26, 28, 31, 40]. However, these lower bounds arc valid asymptotically only 
and hence these methods do not guarantee that the computed approximation is really below 
the exact value. Recently, sufficient conditions for producing lower bounds for eigenvalues 
of symmetric elliptic operators by nonconforming methods have been provided in [19]. The 
described technique is based on satisfying the saturation assumption and on the condition 
saying that the local approximation property of the underlying finite element space have 
to be better than its global continuity property. It is proved that the second condition is 
met by most commonly used nonconforming methods such as the Wilson element, linear 
nonconforming Crouzeix-Raviart element, the enriched nonconforming rotated Qi element, 
the Morley element, the Adini element, and the Morley-Wang-Xu element. The saturation 
assumption is proved for the Morley-Wang-Xu element, the enriched Crouzeix-Raviart ele- 
ment, and the Wilson element. Furthermore, new nonconforming methods satisfying these 
properties are proposed. However, no numerical experiments are presented. 

Further, let us point out a recent result [27], where two-sided a priori bounds for the 
discretization error of eigenvectors are given. 

The method of a priori-a posteriori inequalities that can be used for computation of lower 
bounds on eigenvalues was described and published in [22, 23, 34]. However, in these original 
publications C^-continuous test and trial functions has been used in order to compute the 
actual lower bound. These functions are difficult to work with and therefore, we couple 
the original method of a priori-a posteriori inequalities with the complementarity technique, 
where a certain flux function has to be reconstructed, see, e.g., [2, 14, 18, 32, 36, 38]. This 
flux reconstruction can be done in many ways. We choose the most straightforward approach 
that can be handled by standard Raviart-Thomas flnite element method. 

Moreover, we generalize the original method of a priori-a posteriori inequalities to the 
case of a compact operator between a pair of Hilbert spaces. This generalization is especially 
useful for computing two-sided bounds of the optimal constant in the Friedrichs', Poincare, 
trace, and similar inequalities. It is based on the fact that the optimal constant in these 
inequalities is inversely proportional to the square root of the smallest eigenvalue of the 
corresponding symmetric linear elliptic partial differential operator. 

Further, the generalization we have made enables to set up an abstract framework in the 
Hilbert space setting. The abstract results can be then easily applied to symmetric linear 
elliptic partial differential operators and consequently to the optimal constant in the inequal- 
ities of Fricdrichs'-Poincare type. Furthermore, as a byproduct of the abstract setting, we 
obtain a simple proof of the validity of an abstract inequality of Friedrichs '-Poincare type in 
the Hilbert space setting. The particular choices of the pair of Hilbert spaces, corresponding 
scalar products, and the compact operator then naturally yield the validity of the Friedrichs', 
Poincare, trace, Korn's and other inequalities. 

The main motivation for our interest in two-sided bounds of the constants in the Friedrichs', 
Poincare, trace, and similar inequalities stems from the need of these bounds in a posteriori 
error estimation for numerical solutions of partial differential equations. In particular, the 
existing guaranteed upper bounds on the energy norm of the error utilize a kind of comple- 
mentarity technique, see, e.g., [2, 12, 20, 32, 36, 38]. Estimates of this kind contain constants 
from the Friedrichs', Poincare, trace, and similar inequalities. Optimal values of these con- 
stants are often unknown and therefore suitable approximations have to be used in the error 
estimates. These approximations have to provide upper bounds on these constants in order 
to guarantee that the total error estimator is an upper bound on the error. Moreover, they 
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have to be accurate due to the accuracy and efficiency of the error estimates. 

The method presented in this paper provides accurate upper bounds on these constants. 
In addition, this method naturaUy considers the dependence of the optimal constants on the 
equation coefficients and on the boundary conditions. This dependence can be strong [37] 
and its capturing might be crucial for accuracy and robustness of a posteriori error bounds. 

The rest of the paper is organized as follows. Section 2 introduces a general variational 
eigenvalue problem in the Hilbcrt space setting. It uses the spectral theory of compact op- 
erators to prove several properties of this eigenvalue problem including the existence of the 
principal eigenvalue. In Section 3 we naturally prove the abstract inequality of Friedrichs'- 
Poincare type and show the relation between the optimal constant and the principal eigen- 
value. Further, we briefly describe the Galerkin method that yields an upper bound on the 
principal eigenvalue and concentrate on the method of a priori-a posteriori inequalities and 
on an abstract complementarity result leading to a lower bound on the principal eigenvalue. 
Sections 4-6 apply the abstract results to the case of Friedrichs', Poincare, and trace inequal- 
ity and fully computable two-sided bounds on the optimal constants in these inequalities are 
obtained. Presented numerical experiments illustrate accuracy of the method and depen- 
dence of the optimal constants on a nonhomogeneous diffusion parameter. Finally, Section 7 
draws the conclusions. 



2 Variational eigenvalue problem in the Hilbert space 
setting 

Let V and H be two real Hilbert spaces with scalar products (•, •)y and {■,-)h; respectively. 
The norms induced by these scalar products arc denoted by |j ■ ||y and || • \\h- Further, let 
7 : y — > be a continuous, linear, and compact operator. The center of our interest is the 
following eigenvalue problem. Find G M, e V, Ui ^ such that 

{ui,v)v = \i{^Ui,^v)u Mv^V. (1) 

First, let us show that eigenvalues \i of (1) are positive. 

Lemma 1. Ij Ui Cz V is an eigenvector corresponding to an eigenvalue Xi of (1) then ^ui ^ 
and Ai > 0. 

Proof. Since Ui ^ 0, we have by (1) that ^ \\ui\\y = XiW^fUiW^j. Thus, "fUi ^ and Ai has 
to be positive. □ 

Below we show that eigenvalues Xi and eigenvectors Ui of (1) correspond to eigenvalues 
and eigenvectors of a compact operator, respectively. Consequently, these eigenvalues form a 
countable sequence that can be ordered as Ai < A2 < ■ ■ • . To show this correspondence, we 
deffire a solution operator S : H ^ V . If / e i? is arbitrary then the mapping v 1— )■ (/, ^v)h is 
a continuous linear form on V and, hence, the Riesz representation theorem yields existence 
of a unique element S f G V such that 

{Sf,v)v^{fnv)H yv^v. (2) 

Consequently, the solution operator S is linear and continuous. 

The composition of operators S and 7 is a linear, continuous, and compact operator 
57 : V ^ V , see [33, Theorem 4.18 (f)]. In addition, this operator is selfadjoint, because 
definition (2) yields 



{S^u,v)v = {-fu,jv)H = {iv,ju)h = {S-fv,u)v = {u,Sjv)v 
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Therefore, we can use the Hilbert-Schmidt spectral theorem for 5*7, see [17, Theorem 4, 
Chapter II, Section 3] and obtain that V can be decomposed into a direct sum of two 
subspaces 

V = MekeriSj), (3) 

where ker(S'7) ~ {v G V : Sjv — 0} is the kernel of 6*7 and A4 is generated by all eigenvectors 
of the operator S'7 corresponding to nonzero eigenvalues. Let us recall that Ui G V, Ui ^ 
is an eigenvector of S'7 corresponding to an eigenvalue /ij g R if 

577/4 = ^iiUi. (4) 

Furthermore, the Hilbert-Schmidt spectral theorem implies that the system ui,U2,... of 
eigenvectors corresponding to nonzero eigenvalues of (4) is countable and orthogonal in V . 
A simple consequence of (2) and (4) is the orthogonality of images 71/^, z = 1, 2, . . . \a H. In 
this paper we consider the following normalization of these eigenvectors: 

{■ju^,jUj)h = Vi,j = 1,2, . . . , (5) 

where Sij stands for the Kronecker's delta. 

Now, let us observe that eigenproblems (1) and (4) correspond to each other and, there- 
fore, the spectral properties of the compact operator 57 translate to the properties of the 
eigenproblem (1). 

Lemma 2. Considering the above setting, the following statements hold true. 

1. Number Ai G R is an eigenvalue corresponding to the eigenvector Ui € V of (1) if and 
only if Hi = 1/ Xi is a nonzero eigenvalue corresponding to the eigenvector Ui of the 
operator S'7, see (4)- 

2. The number of eigenvalues Xi of (1) such that Xi < M is finite for any M > 0. 

3. The value Ai = inf ||w||y/||7w|||/ is the smallest eigenvalue of (1). 

Proof. 1. Definition (2) yields identity {'^Ui,^v)H — {S^Ui,v)v for all v G V. Hence, the 
equality (1) can be rewritten as {ui^v)v — Xi{SjUi,v)v, which is equivalent to (4) with 
fii = 1/Xi provided Xi ^ and fii ^ 0. Since Lemma 1 guarantees Ai > for all i 1, 2, . . . , 
the only condition is fii ^ 0. 

2. If we denote the spectrum of S7 by cr(S7) then the compactness of S7 implies that the 
set [£,00) n cr(S7) is finite for any e > 0, see [33, Theorem 4.24 (b)]. The statement follows 
immediately from the fact that Ai = l/fii for fii 0. 

3. Since S'7 is selfadjoint, the Courant -Fischer- Weyl min-max principle implies that 

Hi = sup{(S7U,u)y : \\v\\v = 1} = sup — 7--^ — = sup ,, ..^ 
is finite and it is the largest eigenvalue of the operator S'7. Consequently, 

X -1 ( h^WlY' ■ r Ml 

Ai = /ii = sup II ,1^ = ult -pp (6) 

\v(iV,v^Q Ml ) vev,v^o \\-fv\\% 

is the smallest eigenvalue of problem (1). □ 
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3 Abstract inequality of Friedrichs'— Poincare type 

3.1 The proof of the abstract inequahty 

Properties of eigenproblem (1) ean be utilized in a simple way to derive an abstract inequal- 
ity of Friedrichs'-Poincare type. The Hilbert space versions of the particular Friedrichs', 
Poincare, trace, Korn's and similar inequalities easily follow from this abstract result. For 
examples see Sections 4-6. 

Theorem 3 (Abstract inequality). Let j : V ^ H be a continuous, linear, and compact 
operator between Hilbert spaces V and H. Let Ai be the smallest eigenvalue of problem (1). 
Then 

\\iv\\h < C-,\\v\\v yv&V (7) 

— 1/2 

with C-y = Aj . Moreover, this constant is optimal in the sense that it is the smallest 
possible constant such that (7) holds for all v € V . 

Proof. The validity of the abstract inequality follows immediately from (6): 

hv\\H<K'\Ml v«ey. 

— 1 /2 . 

This inequality holds as equality for v ^ ui and thus, the constant = is optimal. □ 

3.2 Upper bound on the smallest eigenvalue 

The upper bound on Ai can be computed by the standard Galerkin method, which is both 
accurate and efficient [5, 6, 11]. The Galerkin method is based on the projection of the 
eigenproblem (1) into a finite dimensional subspacc C V. We seek eigenvectors S V'^, 
7^ 0, and eigenvalues \^ such that 

{ulv'^)v = \'i{lu1nv'^)H yv'^eV\ (8) 

Considering a basis lySj, j = 1, 2, . . . , in V'^ , we can formulate problem (8) equivalently as 
a generalized eigenvalue problem 

Ay, = X'iMy, 

for matrices A and M with entries 

Ajk = {'^k, (fjW and Mjk = (7<Pfc, 7<Pi)ff , j, k = 1,2, . . . , N. 

The eigenvectors yi S and u,-' e are linked by the relation w'' ~ '^^=iyi-j'Pj ■ The 
generalized matrix eigenvalue problem can be solved by efficient methods of numerical linear 
algebra [7]. 

The Galerkin method for eigenvalue problems is very well understood. The convergence 
and the speed of convergence of this method is established for example in [5, 6, 11]. It is well 
known [11] that the Galerkin method approximates the exact eigenvalues from above, hence 

A, <A^ V^ = l,2,.... 

In particular, the upper bound on the smallest eigenvalue Ai and the corresponding lower 
bound on the optimal constant read 

Ai < Aj and (AJ)-^/^ < C^. (9) 
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3.3 Lower bound on the smallest eigenvalue 

In this part wc concentrate on a computable lower bound on the smallest eigenvalue Ai. 
First, wc formulate an auxiliary result. It states that the images 7^;, i = 1,2,... of the 
orthogonal system of eigenvectors Ui satisfy the Parseval's identity. 

Lemma 4. Let Ui, i ~ 1,2, be the above specified orthogonal system of eigenvectors of 
the operator Sj corresponding to nonzero eigenvalues. Let these eigenvectors be normalized 
as in (5). Finally, let €z V be arbitrary. Then 



i=l 



Proof. Due to the decomposition (3), there exist unique components u;^ e Ai and u° G 

u°. Since Sju° = 0, we have = {S^u°,u°)v = hu°,ju°)H 



,M 



,M 



ker(5'7) such that 

by (2) and hence 71/° = 0. Consequently, 7U* = 7U; 

System "fUi, i = 1,2, . . . , forms an orthonormal basis in A4. Thus, we can use the standard 
Parseval's identity in Ai [41, Theorem 2, Chapter III, Section 4] to obtain 



i=l 



□ 



The derivation of the lower bound on Ai is based on the method of a priori-a posteriori 
inequalities. This method relies on an abstract theorem proved in [22]. Wc formulate this 
theorem in the setting of this paper and for the readers' convenience we present its brief 
proof. Notice that in contrast to [22], Theorem 5 operates with a pair of Hilbert spaces and 
with a compact operator between them. 

Theorem 5. Let ^ : V ^ H be a continuous, linear, and compact operator between Hilbert 
spaces V and H . Let £ V and A* e M 6e arbitrary. Let us consider w &V such that 



If 7W* 7^ then 



{w,v)v = iu*,v)v - X*i'ju*,^v)H Vw e K 

\iw\\h 



A,: - A, 



A. 



< 



\nu*\\H 



(10) 

(11) 



Proof. Using Lemma 4 and definitions (1) and (10), we obtain 



A,; — Ai 



l7"*llff <E 

00 

= E 

00 

= E 



A,; — A, 



A,; 



{w,Ui)v 



A,: 



A,; 



E \{lW,-iUi)H\^ = WlwfH- 



□ 
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In order to obtain a computable lower bound on Ai, we combine the estimate (11) with 
a complementarity technique, see Sections 4-6. The bounds derived by the complementarity 
technique depend on the particular choice of spaces V and H, but they have a common 
general structure. The following theorem utilizes this general structure and presents an 
abstract lower bound on Ai. 

Theorem 6 (Abstract complementarity estimate). Let G A* G K be arbitrary and let 
w € V satisfy (10). Let Ai be the smallest eigenvalue of (1) and let the relatively closest 
eigenvalue to A* be Xi, i.e. let 



Ai — A* 


< 


A; — A* 


Ai 




A, 



Vi = l,2,.... (12) 

Finally, let us consider A> and B > such that 

\\w\\v < A + C^B, (13) 
where B < X*\\"/u^\\h and Cj is the optimal constant from (7). Then 

xi < Ai and < I/X2, (14) 



where 



^(-«+V«' + 4(A, -/?)), a = — and p = (15) 



Proof. Using the fact that = Aj^ , the estimate (11), assumptions of the theorem, and 
inequality (7), wc obtain the validity of the following relation 



^* - ^1 ^ ______ Aj - A, ||7w||g \\w\\v ^ r< ^^r^R 



A*C — 1 = — < min 

Ai i 



This is equivalent to the quadratic inequality 



0<C2(/3-A,) + C^a + l. 

Solving it for C-y under the assumption < A,, wc conclude that this inequality is not 
satisfied for > I/X2. Thus, has to be at most I/X2. □ 

The particular complementarity estimates have the form (13), where the numbers A and 
B arc obtained by an approximate minimization procedure, see Sections 4-6. The better ap- 
proximation of the exact minimizer we compute the tighter bound (13) and consequently (14) 
we obtain. The exact minimizer yields equality in (13) and B = 0. Therefore, the assumption 
B < A*||7u*j|H can always be satisfied by computing sufficiently accurate minimizer. 

The assumption (12) is crucial and it cannot be guaranteed unless lower bounds on Ai 
and A2 arc known. However, since the Galerkin method is known to converge [5, 6, 11] 
with a known speed, very accurate approximations of Ai and A2 can be computed. If these 
approximations are well separated then they can be used in (12) to verify its validity with a 
good confidence. 

In order to increase this confidence, we propose a test. This test is based on the following 
observation. Let A'i°" < Ai < A* < X^^"" < A2 < A2'', Di = (A* - A'i°*)/A'i°" and D2 = 
{X^2™ ^ A*)/A2''. Then inequality Di < D2 implies the assumption (12). Indeed, if all 
these inequalities are satisfied then, clearly, |(Ai — A*)/Ai| < Di < D2 < |(A2 — A*)/A2| < 
|(Ai - A*)/A,| for alH = 2, 3, . . . and (12) holds true. 
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Thus, if we knew guaranteed lower bounds ^-i^d A2 on the first eigenvalue Ai and 
the second eigenvalue A2 of (1), respectively, then we could use the Galcrkin method to 
compute upper bounds A* and X^^ , check if A* < Aj"™ and the inequality Di < D2 would 
then guarantee the validity of (12). However, the guaranteed lower bounds A'™ and Aj"^ are 
not available. 

Therefore, in practice we propose to set X^^^ = X|, see (14), A^""^ = (A, + A2^)/2, and 
compute Di and D2 with these values. This yields the following diagnostics indicating the 
validity of assumption (12). li D2 < ot Di > D2 then it is highly probable that some of 
the assumptions is not satisfied and the results should not be trusted. On the other hand 
if D2 > and Di is several times smaller than D2 then we can have a good confidence 
about the validity of (12). This diagnostics performs very well in all numerical experiments 
presented below. At the early stages of the computations the approximations are not very 
precise and the diagnostics showed that the results are untrustworthy. However, at the final 
stages of the computations the value Di was at least ten times smaller than D2 providing 
good confidence about the validity of (12). 

4 Application to the Friedrichs' inequality 

In this section we apply the above general theory to the case of Friedrichs' inequality. We will 
consider the variant of Friedrichs' inequality that is suitable for general symmetric second- 
order linear elliptic differential operators. First, we introduce differential operators of this 
general type in Subsection 4.1. In Subsection 4.2 we derive two-sided bounds on the opti- 
mal constant in the Friedrichs' inequality that are based on the general theory and on the 
complementarity technique. In Subsection 4.3 certain computational issues are discussed and 
finally, in Subsection 4.4 we present numerical results. 

4.1 A general symmetric second-order linear elliptic operator 

Let us consider a domain i7 C M"* with Lipschitz boundary 9f2. Let dil consist of two 
relatively open parts Fd and Fn such that dD. = Fd U Fn and Fd n Fn =0- Note that we 
admit the case where either Fd or Fn is empty. Further, let us consider a matrix function 
A e [L°°{n)]'^^'^, coefficients c G L°°{fl), and a e L°°(Fn). Matrix A is assumed to be 
symmetric and uniformly positive definite, i.e. there exists a constant C > such that 

i^Aix)$, > Clip V| e R'' and for a.e. x e n, 

where | • | stands for the Euclidean norm. Coefficients c and a arc considered to be nonnegative. 
Further, we consider a subspace 

i/pj^(51) = {v G H^{il) ; w = on Fd in the sense of traces} 

of the Sobolev space H^{n) of L'^{n) functions with square integrable distributional deriva- 
tives. In what follows, we use the notation (•,•) and (^^^^(rN) fo'^ the L^(fi) and ^^(Fn) 
scalar products, respectively. Using this notation, we define a bilinear form 

a{u,v) = {AWu,Wv) + {cu,v) + (au,-(;)i2(rN)- (16) 

This bilinear form is a scalar product in H^_^{il) under the conditions presented in the 
following lemma. Its proof follows for instance from [21, Theorem 5.11.2]. 

Lemma 7. The bilinear form defined in (16) is a scalar product in H^^{il) provided that at 
least one of the following conditions is satisfied: 
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(a) mcaSd„irD > 0, 

(b) there exists a nonempty ball B d D, such that c > on B, 

(c) there exists a subset Ta C Fn such that mcasd-i > and a > on F„. 

The equivalence of the norm induced by a(-,-) and the standard norm yields com- 
pleteness and therefore the space Hp^{i}) equipped with the scalar product a(-, •) is a Hilbert 
space. This enables to use the theory from Section 2. We set 

V^H^jn), {u,v)v ^a{u,v), H = L\n), {u,v)h = M, (17) 

and we define 7 : H^^{il) — > L^(r2) as the identity, which is compact due to the Rellich 
theorem [1, Theorem 6.3]. With this setting and with the notation || • ||a for the norm 
induced by a(-, •), we obtain the validity of the Friedrichs' inequality. 

Theorem 8. Let the bilinear form a(-, •) given by (16) form a scalar product in Hp [il). 
Then there exists a constant Cp > such that 

\\v\\mn) < CF\\v\\a yvGH^jn). (is) 

— 1/2 

Moreover, the optimal value of this constant is Cp = A-^ , where Ai is the smallest eigen- 
value of the following problem: find Ui € Hp^{rt), Ui 7^ 0, and A,; G M such that 

aiu„v) = X,{u,,v) yveH^jn). (19) 

Proof. Lemma 7 guarantees that (51) is a Hilbert space with scalar product a(-, •) given 
by (16). The statement then follows immediately from Theorem 3. □ 

Let us note that the most common version of Friedrichs' inequality 

\\v\\L-(n) < CF\\Vv\\L2^n) yveH^in) 

follows from (18) with Fd = 957, Fn = 0,-4 being identity matrix, c = 0, and a = 0. As 
usual, we denote by 11^(^1) the space iJp^(r2) with Fd = dft. 

4.2 Two-sided bounds on Friedrichs' constant 

A lower bound on the Friedrichs' constant Cp can be efficiently computed by the Galerkin 
method. We use the setting (17) and proceed as it is described in Section 3.2. The up- 
per bound on Cp is obtained by the complementarity technique presented in the following 
theorem. 

Let H(div, O) stands for the space of d-dimensional vector fields with square integrable 
divergence and let n be the unit outward-facing normal vector to the boundary dfl. Further, 
let ||q||3i = {Aq,q) be a norm in [L'^{il)]'^ induced by the matrix A. 

Theorem 9. Let V = (51), u* G V and A, G K. Let the bilinear form a{-, •) given by 
(16) form a scalar product in Hy^{VL). Let w £V satisfy 

a(w, w) = a(it,, w) — A,(u*, w) Vw G T^, (20) 

where recall that (•, •) stands for the L^(51) scalar product. Then 

||«^i|a < ||Vu* - + Cp||A,w, - cu, + divq||i2(n) Vq G W", (21) 

where W = {q G H(div, f2) : q • n = — au* on Fn}- 
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Proof. Let us fix any q G W, test (20) hy v = w and use the divergence theorem to express 

\\w\\l = {AVu*,Vw) + {cu^.,w) + (au*, u;)L2(rjj) - A*(u*,ii;) - {q,Vw) - {A\v q.w) 

+ {q -n, w)i2(rN) = (^(Vu* - A'^q), Vw) - (A*m* - cu* + divq, li)). 

The Cauchy-Schwarz inequahty and Friedrichs' inequahty (18) yield 

\\w\\l < 11^"* - •^^^q\\a\\'^w\\a + CpjlA^u* - cu* + div(7||i2(s2)|lw|la. 
Inequality ||Vi/7j|^ < \\w\\a finishes the proof. □ 
The estimate (21) is of type (13) with 

WWu^ - A^^q\\A, B = ||A,M* - cu* + divq||i2(n), (22) 

and C-y = Cp. The numbers A and B can be readily computed as soon as suitable approx- 
imations A*, M*, and a suitable vector field q £ W arc in hand. Estimate (14) then gives a 
guaranteed upper bound on the Friedrichs' constant Cp. 

The crucial part is the computation of suitable approximations A, and of the smallest 
eigenvalue Ai of problem (19) and its corresponding eigenvector ui such that the inequality 
(12) is satisfied. Equally crucial is a suitable choice of the vector field q G W in such a way 
that estimate (21) provides a tight upper bound on lluijla- A possible approach and practical 
details about these issues are provided in the next section. 

4.3 Flux reconstruction 

In order to compute two-sided bounds on Friedrichs' constant, we proceed as follows. First, 
we use the Galerkin method, see Section 3.2, to compute an approximation A^ and Ui of 
the eigenpair Ai and ui of (19). We set A* = X'l, m, = it^, compute a suitable vector field 
q £ W, and evaluate the numbers A and B by (22). Estimate (14) then provides the upper 
bound on the Friedrichs' constant Cp and the Galerkin approximation A^* yields the lower 
bound (9) on Cp. 

The key point is the computation of the suitable vector field q. For simplicity, we choose a 
straightforward approach of approximate minimization of the upper bound (21) with respect 
to a suitable subset of W. First, we exploit the affinc structure of W. Let us choose an 
arbitrary but fixed q G W. The practical construction of this q is a geometrical issue 
depending on 51 and coefficient a. It suffices to construct a vector field qi G H(div,17) such 
that • n = 1 on Fn and a function a G H^(Q,) such that a = a on Fn. Let us note 
that usually Raviart-Thomas-Nedelec space is considered for construction of vector fields in 
H(div, ri). Then we can simply set q = —au^qi. This q obviously satisfies the boundary 
condition q ■ n = — au* on Fn and if qi and a are sufficiently smooth then q G H(div, 57). 

In any case, having a q G W^, we can express the affine space W as W = q + Wo, where 

Wo = {q G H(div, r2) : q • n = on Fn} 

is already a linear space. The idea is to minimize the upper bound (21) over the set q + Wq, 
where Wq C Wq is a finite dimensional subspace. However, the right-hand side of (21) is a 
nonlinear functional in q and, thus, in order to simplify the computation, we use the idea 
from [32] and approximate it by a quadratic functional. We rewrite inequality (21) using the 
notation (22) and apply an elementary estimate: 



Ml <{A + CpB f < (1 + Q-^)A^ + + q)CIB^ > 0. 
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The right-hand side of this inequahty is aheady a quadratic functional for a fixed g, but the 
exact value of Cp is unknown in general. However, it is sufficient to find an approximate 
minimizer only. Therefore we approximate Cp by the available value (A^*)"^/^. This leads us 
to the minimization of 

(1 + g-')\\Vu1 A-^qfA + (1 + g){X'D-'\\^i^i cu^i + divglli.(a) (23) 

over the affine set q+W^ with a fixed value oi g > 0. This minimization problem is equivalent 
to seeking Qq G Wq satisfying 

Biq'^,w'^)^Tiw'^)-Biq,w'$) yw'$eWo\ (24) 

where 

B{q,w) = (divq, divu;) + —{A^^q.w), 

e 

T{w) = -±{\Jv!l,w) - (A>^ -c4,div«;). 
Q 

The computed vector field q'' = q + e is then used in (22) to evaluate A and B and 
consequently the two-sided bounds C|?™ < C-p < Cp^, where 

= (Aj)-i/2 and CpP = l/X2, (25) 

X2 is given by (15) with = , see (9) and (14). 

Note that problem (24) can be naturally approached by standard Raviart-Thomas finite 
elements [13]. Further we note that a particular value of the constant g can influence the 
accuracy of the final bound. However, this infiuence was minor in all cases we numerically 
tested and the natural value g — 1 yielded accurate results. If necessary a simultaneous 
minimization of (23) with respect to both g> Q and q Gq + Wq can be performed. 

4.4 Numerical experiment 

In order to illustrate the capabilities of the above described approach for computation of two- 
sided bounds on Friedrichs' constant, we present numerical results showing the dependence of 
Friedrichs' constant on piecewise constant values of A. We consider the general setting from 
Section 4.1 and in particular we set = (—1, 1)^, Fn = {{xi,X2) G : xi = 1 and — 1 < 
X2 < 1}, Fd = dfl \ Fn, c 0, and a = 0. The matrix A is piecewise constant, defined as 
A{xi,X2) = I for X1X2 < and A{xi,X2) = al for X1X2 > 0, where / € R'^^'' stands for the 
identity matrix and the value of the constant a is specified below. 

We employ the standard lowest-order (piecewise linear) triangular finite elements to dis- 
cretize eigenvalue problem (19). Thus, we consider a triangular mesh and seek Galerkin 
solution of problem (19) in space 

V'' = G H^jn) : v'^Ik e P\K), VA- e %}, (26) 

where P^{K) denotes the space of affine functions on triangle K, see Section 3.2. We point 
out that the discontinuity of A{xi,X2) causes a strong singularity of the eigenvectors at the 
origin. Therefore, we employ a standard adaptive algorithm, see Algorithm 1, and construct 
adaptively refined meshes in order to approximate the singularity well. We use the localized 
version of (23) to define the error indicators 

V% = (1 + - -A-'qWXk + (1 + g){X'l)-'\\X'M cul + divq||2,(^^ (27) 
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Adaptive algorithm 
Step 1 Construct an initial mesh Th- 

Step 2 Use the space (26) and compute Galerkin approximations Aj e M and u'^ G V'^ of 
the smallest eigenvalue and the corresponding eigenvector of problem (19), sec (8). 

Step 3 Use the space of second-order Raviart-Thomas [13] finite elements 
= {wh e Wo : Wh e [P2(K)]^ V/^ e %} and solve (24). 

Step 4 Evaluate two-sided bounds (25). Set C^^ = (Cp^ + C^™)/2, 

£^REL = (Cp^ ~ C'f*)/Cp^^, and stop the algorithm as soon as Srel < Etol- 

Step 5 Compute error indicators (27) for all K ^ Th and sort them in descending order: 
VKi ^ '7^2 ^ ■ ■ ■ ^ '?A']v [ I where A^oi is the number of elements in Th ■ 

Step 6 (Bulk criterion.) Find the smallest n such that 6'^ Y^f=i Vk^ ^ '7i'ii where 

6* e (0, 1) is a parameter. 

Step 7 Construct a new mesh Th by refining elements Ki, K2, . . . , Kn- 

Step 8 Go to Step 2. 

Algorithm 1: Mesh adaptation algorithm for two-sided bounds on Friedrichs' constant. 



d 




L-F 


£^REL 


^DOF 


0.001 


9.0086 


9.0939 


0.0094 


4 832 


0.01 


2.8697 


2.8971 


0.0095 


5 003 


0.1 


1.0035 


1.0124 


0.0088 


7866 


1 


0.5693 


0.5743 


0.0086 


4 802 


10 


0.3173 


0.3201 


0.0088 


7866 


100 


0.2870 


0.2897 


0.0095 


5 003 


1000 


0.2849 


0.2876 


0.0094 


4 832 



Table 1: Friedrichs' constant. The lower bound Cp^ , upper bound Cp^, relative error -Erel, 
and the number of degrees of freedom iVDOF = dim for particular values of d. 

for all K G %, where = {■Aw,w)i^2(^j^y 

The unknown value of Friedrichs' constant lies between bounds Cp™ and Cp^ computed 
in Step 4 of Algorithm 1. These bounds are utilized for the relative error stopping criterion. 
When the algorithm stops, the relative error is guaranteed to be at most the given tolerance 
-E'tol- 

In this particular numerical experiment we have chosen the initial mesh with eight tri- 
angles as shown in Figure 1. The marking parameter in Step 6 and the tolerance for the 
relative error in Step 4 were chosen as 6 ^ 0.75 and -Etol = 0.01. Further, we naturally set 

The results for a series of values of d are summarized in Table 1. For each particular value 
of a, we run the adaptive algorithm until the relative error drops below Etol = 0.01. This 
error level was reached in all cases using several thousands of degrees of freedom. Notice the 
considerable dependence of the optimal value of Friedrichs' constant on d. The values for 
d = 0.001 and d = 1000 differ more than thirty times. Further notice that the exact value of 
Friedrichs' constant for a = 1 is Cp = A/{tt\/5) ~ 0.5694. 
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(-1,1) Td (1,1) 




(-1,-1) Td (1,-1) 

Figure 1: The initial mesh for the adaptive algorithm. The Dirichlet and Neumann parts 
of the boundary are indicated as well as the piecewise constant matrix A: A = I in white 
elements and A = al in gray elements. 



5 Application to the Poincare inequality 
5.1 Poincare inequality and the optimal constant 

In this section we consider the case when none of conditions (a)-(c) of Lemma 7 is satisfied. 
Therefore, we assume the general symmetric second-order elliptic operator as described in 
Section 4.1 with c = 0, a = 0, Fd = 0, and Fn = dil. We apply the general theory of 
Section 3 with V = H^in) = {w e H\n) : {v,l) 0}, {u,v)v = a{u,v) = (^Vw,Vw), 
H = L'^{n), and {u,v)h — iu,v). It is an easy exercise to verify that a(-,-) forms a scalar 
product on H^{rt) and that it induces a norm || • ||a equivalent to the standard iJ^-seminorm 
and _ff^-norm. The operator 'y : V H is simply the identity mapping / : H^{H) — > 
This mapping is clearly compact, because the identity mapping from if ^($7) to H^{il) is 
linear and continuous and the identity mapping from H^(fl) to L^(fl) is compact due to the 
Rellich theorem [1, Theorem 6.3]. This setting enables to use the general result of Theorem 3 
and obtain the following result. 

Theorem 10. There exists a constant Cp > such that 

\\v\\L-m <Cp\\v\\a yveH\n). (28) 

— 1/2 

Moreover, the optimal value of this constant is Cp = Aj , where Ai is the smallest positive 
eigenvalue of the following problem: find Ui G H^{n), Ui ^ 0, and A^ € M such that 

a{u„ v) = A,(m„ v) Vw e H\n). (29) 

Proof. The inequality (28) follows immediately from Theorem 3. This theorem also implies 

~ —1/2 

that the optimal constant is Cp = Aj^ ' , where Ai is the smallest eigenvalue of the following 
problem: find Ui G H^{il), Ui ^ 0, and Ai G K such that 

a{u,,v) = \{u„v) yveH^n). (30) 

Notice that < Ai < A2 < • ■ • • Similarly, the eigenvalues of (29) satisfy = Aq < Ai < A2 < 
. . . and the zero eigenvalue corresponds to a constant eigenvector uq = 1. It can be easily 
shown that Ai = A^ and Ui = Ui for alH = 1, 2, . . . . Thus, the smallest eigenvalue Ai of (30) 
is equal to the smallest positive eigenvalue Ai of (29) and the proof is finished. □ 
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Let us note that the Poincare inequaUty (28) can be equivalently formulated as 

\\v~v\\L'^in)<Cp\\v\\a yveHHn), i^{v,i)/\n\. (31) 

Further, it is clear that the common version of the Poincare inequality 

\\v - vh-in) < Cp\\Vv\\L2^n) G H^{n) 
follows from (31) and consequently from (28) with A being the identity matrix. 

5.2 Two-sided bounds on the Poincare constant 

The lower bound on the Poincare constant Cp can be computed in the standard way by 
the Galerkin method, see Section 3.2. The only difference is that here we compute the 
approximation Xi of the second smallest (the smallest positive) eigenvalue of (29), because 
the smallest eigenvalue is Aq ~ 0. In order to compute the upper bound, we employ the 
complementarity technique as follows. 

Theorem 11. Let V = H^{^}), e V, and A, e R. Let w eV satisfy 

a{w,v) ~ a{u^.,v) ~ X^.{u^,,v) Vv e V. (32) 

Then 

i|w;||a < ||Vu, -^"iq|U + Cp||A*u* +divq||i2(s^) Vq e Wq, (33) 
where Wq = {q G H(div, fi) : q • n = on dft}. 

Proof. Fixing arbitrary q e Wq, testing (32) hy v = w and using the divergence theorem, we 
obtain 

= (^Vit*, Vw) - A*(u*,w) - {q,Vw) - (divq,u)) 

= (^(Vu* - A^^q),Vw) - (A*'u* + divq, w). 
The Cauchy-Schwarz inequality and Poincare inequality (28) yield 

\\w\\l < ||Vm* -^"^q|U||Vw|U + Cp||A*M* + divq||i2(o)||u>i|a. 
Since ||Vwj|^ = the proof is finished. □ 

We observe that complementarity estimate (33) is of type (13) with 

\\\/u^ - A^^q\\A, S = ||A,u, +divq||i2(o), (34) 

and = Cp. As soon as a suitable vector field q € Wq is available, the numbers A and B 
can be computed and used in (14)-(15) to obtain a guaranteed upper bound on the Poincare 
constant Cp. 

A straightforward approach for computing a suitable vector field q'* S Wq was described in 
Section 4.3 for the case of Friedrichs' constant. It can be directly used also for the Poincare 
constant. It is even simpler, because c = 0, a = 0, and Pn = dfl. The only difference 
is that the approximation A^* of the second smallest (the first positive) eigenvalue of (29) 
has to be used. In particular, the approximation q'' e Wq is computed using (24). This 
vector field is then used in (34) to evaluate A and B and consequently the two-sided bounds 
C|?^ < Cp < CpP, where 

= (Af)-i/2 and CpP = l/X2, (35) 

X2 is given by (15) with 711* = Ui being the approximate eigenvector of (29) corresponding 
to Aj-, see also (9) and (14). 
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a 






Erel 


^DOF 


0.001 


14.2390 


14.3690 


0.0091 


3400 


0.01 


4.5199 


4.5623 


0.0093 


3510 


0.1 


1.4849 


1.4989 


0.0094 


4 382 


1 


0.6365 


0.6424 


0.0092 


3 009 


10 


0.4696 


0.4740 


0.0094 


4 382 


100 


0.4520 


0.4562 


0.0093 


3510 


1000 


0.4503 


0.4544 


0.0091 


3400 



Table 2: Poincare constant. The lower bound Cp™, upper bound Cp^, relative error i^REL, 
and the number of degrees of freedom A^dof = dim V'^ for particular values of a. 

5.3 Numerical experiment 

We consider the same setting as in Section 4.4. The only difference is that in the case of 
Poincare constant we assume ~ dil and = 0. We use the adaptive algorithm as before 
(Algorithm 1) with clear modifications. We use error indicators (27), where A^* and are 
the Galerkin approximations of the second smallest (the smallest positive) eigenvalue Ai of 
(29) and its corresponding eigenvector ui. The relative error in Step 4 of Algorithm 1 is 
computed using two-sided bounds (35). 

The obtained two-sided bounds for a series of values of 5, are presented in Table 2. The 
guaranteed 0.01 relative error tolerance was reached in all cases using several thousands 
degrees of freedom. As in the case of Friedrichs' constant, we observe considerable dependence 
of the Poincare constant Cp on a. Finally, we point out that the exact value of the Poincare 
constant for a = 1 is Cp = 2/7r w 0.6366. 

6 Application to the trace inequality 
6.1 Trace inequality and the optimal constant 

In order to apply the general theory from Sections 2-3 to the case of the trace inequality, 
we consider the same general symmetric second-order elliptic operator as in Section 4.1. In 
addition we assume meas^-i Fn > 0. The general theory is applied with V = H^^{Q), 
{u,v)v = a{u,v), H ~ L^(rN), and {u,v)h = {u,v)]^2(Y^y The operator j : V ^ H is the 
standard trace operator. Its compactness and other properties are provided in [21, Theorem 
6.10.5], see also [10]. The general result from Theorem 3 then translates as follows. 

Theorem 12. Let the bilinear form a(-, •) given by (16) form a scalar product in (fl). 
Then there exists a constant Ct > such that 

MlL^(r«) < CtMU yveH^jn). (36) 

— 1/2 

Moreover, the optimal value of this constant is Ct = , where Ai is the smallest eigen- 

value of problem: find Ui G Hp_^{fl), Ui ^ 0, and A, G M such that 

a{u„v) = Xi{ui,v)L2(^r^-j \/v e Hl^{n). (37) 

Proof. The statement follows immediately from Lemma 7 and Theorem 3. □ 

Clearly, the common version of the trace inequality 

\\v\\man)<CT\\Vv\\mn) e H\n) 

follows from (36) with Fq = 0, Fn = 9ri, A being the identity matrix, c = 0, and a = 0. 
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6.2 Two-sided bounds on the trace constant 

As in the case of Friedrichs' inequality, we eompute the lower bound of the optimal value 
for the constant Ct by the Galerkin method, see Section 3.2. In order to compute an upper 
bound on Ct we employ the complementarity technique as follows. 

Theorem 13. Let V = ffp (fl), u^, G V and A» G M. Let the bilinear form a(-, •) given by 
(16) form a scalar product in V . Let w &V satisfy 

a{'w,v) = a{u^.,v) - K{u^,v)l2(j.^) Mv ^V. (38) 

Then 

\\w\\a < |lV7i, - A~'^q\\A + Cf\\cu^ - divq'|li2(f2) + CT\\au^ - Ku^ + q ■ n|li2(rN) (39) 
for all q e H(div, il). 

Proof. Let us fix any q G H(div, ft), test (38) hy v = w and use the divergence theorem to 
express 

ll«^lla = (-4Vw*, Vw) + {cut,w) + (Q;uH.,u;)L2(rN) " A,(u,, u;)L2(r„) 
- (q, Vw) - (divq,-u;) + (q • n,w)L2{r^) 
= {A{Vur — A^^q),Vw) + {cu^, — div q, w) + (au, — A*u, + q ■ n, w)l2^y^). 

The Cauchy-Schwarz inequality, Friedrichs' inequality (18), and trace inequality (36) yield 

||w|la < l|Vw* - ^"^g|U||Vw|U + CfIIcuh. - divq||o!|w!|a 

+ Cxllau* - A*w* + q ■ n\\L2^Y-^^\\w\\a. 

Inequality ||Vw||^ < ||w||a finishes the proof. □ 

As in the case of Friedrichs' inequality, the bound (39) is of the type (13) with 

A = ||Vit, - A^^q\\A + CpjlCTt, - divqj|i2(jj), B = \\au^ - A*u* + q ■ ni|i2(r„) (40) 

and C-y = Ct. Let us note that the complementarity estimate (39) is just one out of several 
possibilities. This bound comes from [32] and contains the Friedrichs' constant Cf. Instead 
of its exact value an upper bound as computed in Section 4 can be used here. However, there 
exist other variants of the complementarity technique that can be used to obtain a bound on 
\\w\\a without the Friedrichs' constant Cp. See for example [2, 12, 35, 38]. 

A suitable vector field q € H(div, VI) is computed by approximate minimization of the 
right-hand side of (39). In a similar way as we obtained the functional (23), we obtain the 
quadratic functional 

(1 + g-'){l + a-i)|l V4 - A-'q\\\ + (1 + e){l + a-^){Cl^ncu\ divqH^.^^, 

+ {l + a){\1)-^\\au1-\1u1 + q-n\\l2^-Y^) > 0, a > 0, (41) 

where Cp^ is an upper bound of Cf computed as described in Sections 4.2-4.3, Aj is the 
approximation of the smallest eigenvalue of (37) obtained by the Galerkin method and u\ is 
the corresponding approximate eigenvector. We look for the minimum of (41) over a finite 
dimensional subspace of H(div, VI). This minimization problem is equivalent to seeking 

B{q'\w)^F{w) -iweW'', 
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a 






Erel 


^DOF 


0.001 


17.8110 


17.9760 


0.0092 


5 523 


0.01 


5.6490 


5.7047 


0.0098 


5 418 


0.1 


1.8433 


1.8593 


0.0086 


7 775 


1 


0.7963 


0.8033 


0.0088 


5 499 


10 


0.5829 


0.5880 


0.0086 


7 775 


100 


0.5649 


0.5705 


0.0098 


5 421 


1000 


0.5632 


0.5685 


0.0092 


5 523 



Table 3: Trace constant. The lower bound CJ^™, upper bound C^^ , relative error £'rel, and 
the number of degrees of freedom A^dof = dim V'^ for particular values of d. 

where 

Biq,w) = ii^A^(C^Pf(divq,divt^) + l^X'l{A-'q,w) + (q ■ n,w ■ n)i2(r,), 
:F{w) = l^X'l(yu1,w) + ^X'liC^nHculdivw) + (A?4 - aulw n)L2(r,). 

Practically, the classical Raviart-Thomas finite element method [13] can be used to solve 
this problem. Note that the natural values for g and a are g = 1 and cr = 2, because then 
(1 + g^^)il + (T^i) = (1 + g)il + cr^i) = (1 + (t) = 3, see (41). Similarly to the case of 
Friedrichs' constant, these natural values often yield accurate results. If not, a simultaneous 
minimization of (41) with respect to p > 0, cr > 0, and q G can be performed. 

The computed vector field q^ G is then substituted to (40) to get A and B and 
consequently the two-sided bounds < Ct < C^^, where 

C^™ = (A^)-i/2 ^j^j CtP = 1/X2, (42) 
X2 is given by (15) with — u^, see (9) and (14). 



6.3 Numerical experiment 

To illustrate the numerical performance of the above described method, we consider the 
same example as in Section 4.4. We proceed in the same way as in Section 4.4 with clear 
modifications in order to compute two-sided bounds (42) on the trace constant. As before, 
the adaptive algorithm is steered by error indicators that are defined by a localized version 
of (41): 

ril^{l + Q-'){1 + o-'mu\ A-'q\\\j, + {l + g){l + o-^){Cl^ncu\ divq||i.(^,) 

+ (1 + a){\1)-^\\au1 - Xlul + q ■ n||i.(a^,nr.) G Th- 

The parameters g and a are naturally chosen as gi = 1 and a ~ 2. The values of the upper 
bound are taken from Table 1. 

The obtained results are presented in Table 3. The method succeeded in obtaining guar- 
anteed two-sided bounds on the trace constant with a relative error at most 0.01 in all 
cases using several thousands of degrees of freedom. The particular value of the trace con- 
stant Ct depends considerably on a. The values for a = 0.001 and a = 1000 differ more 
than thirty times. Finally, notice that the exact value of the trace constant for a = 1 is 
Ct = (2/(7rcoth7r))i/2 « 0.7964. 



TWO-SIDED BOUNDS FOR EIGENVALUES 



18 



7 Conclusions 

We present a method for computing guaranteed lower and upper bounds of principal eigen- 
values of elliptic operators and consequently for computing guaranteed two-sided bounds of 
the optimal constants in Friedrichs', Poincare, trace, and similar inequalities. The bounds 
are guaranteed provided there are no round-off errors and all integrals are evaluated exactly. 
Further, the bounds are guaranteed only if the domain J7 is represented exactly by used finite 
elements. Furthermore, the upper bounds on eigenvalues computed by the Galerkin method 
are guaranteed only if the corresponding matrix-eigenvalue problems are solved exactly. On 
the other hand, the lower bounds on eigenvalues obtained by the complementarity technique 
are guaranteed even if the matrix-eigenvalue problems and linear algebraic systems are solved 
approximately only. In any case, the crucial assumption for having guaranteed lower bounds 
on eigenvalues is (12). 

These two-sided bounds can be of interest if the corresponding eigenvalue problem cannot 
be solved analytically and if analytical estimates are not available or they are too inaccurate. 
In particular, this is the case of complicated geometry of the domain Jl, mixed boundary 
conditions, presence of non-constant and/or anisotropic diffusion coefficient A, presence of 
reaction coefficient c, and presence of coefficient a. 

The method is quite general and it can be used for a wide variety of problems. The general 
Hilbert space setting presented in Section 2 enables a variety of applications including linear 
elasticity. We believe that this approach can be further generalized. Nonlinear eigenvalue 
problems and nonsymmetric operators can be of particular interest and generalizations in 
these directions can be subject for further research. 
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